Universal relation between longitudinal and transverse 
conductivities in quantum Hall effect 

Igor Ruzin 1 and Shechao Feng 1 ' 2 
1. Department of Physics, Unversity of California, Los Angeles, CA 90024-1547, USA; 
2. Department of Physics, Hong Kong University of Science and Technology, 
Clear Water Bay Road, Kowloon, Hong Kong 
(February 1, 2008) 



^t- : 

On ■ 

s ■ 

<\ 
t> ■ 



> 

O 

oo 
O 

ON 



c 

o 
o 



We show that any critical transition region between two 
adjacent Hall plateaus in either integer or fractional quan- 
tum Hall effect is characterized by a universal "semi-circle" 
relationship between the longitudinal and transverse conduc- 
tivities a X x, o~ xy , provided the sample is homogeneous and 
isotropic on a large scale. This is demonstrated both for the 
phase-coherent quantum transport as well as for the incoher- 
ent transport. 
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Recently there has been much interest in the study 
of transition regions between successive quantum Hall 
plateaus, in both the integer and the fractional regime. 
Most of these studies have been focused on how the tran- 
sition region width AB scales with temperature at low 
temperatures. Experimentally, it has been noted that 
the critical dependence AB oc T K can be extracted ei- 
ther directly from the half- width of the peak in a xx {B), 
or from the sharpness of the step in Hall conductivity 
a xy (B). Thus it has been realized, on intuitive ground 
at least, that the peak- like o~ xx (B) function and the step- 
like a xy (B) function in the transition region between two 
successive plateaus are related. Based on a semiclassical 
model that the transition region, in the critical regime, 
is represented by a random mixture of two liquids with 
different quantized local Hall conductivities o~\ and a 2 , 
Dykhne and Ruzin [Q recently developed a phenomeno- 
logical theory to show that a xx and o~ xy are related by a 
simple "semi-circle" relation: 
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Physically, a finite effective longitudinal conductivity a xx 
at T — > is related to the existence of singularities in the 
current distribution ("hot spots") close to the percola- 
tion threshold. To derive relation ([!]), authors in Ref. 
assumed a very weak scattering in the bulk and described 
it by a local diagonal conductivity a <C cr 2 — ax, defined 
on a scale shorter than the correlation length of the two- 
phase mixture. This assumption is essentially classical 
and does not apply to the more experimentally relevant 



low-temperature regime where quantum coherence exists 
on a much larger scale . 

In this letter, we demonstrate that the same simple 
relationship between a xx and o~ xy holds in the quantum 
coherent transport regime as well. We achieve this by 
first proving a theorem that the semi-circle relation be- 
tween a xx and a xy is equivalent to the perpendicularity 
between the average current densities in the two phases. 
We then show by simulations on the quantum percola- 
tion model that the currents are indeed perpendicular to 
a high accuracy, thus finishing the demonstration. 

As in Ref. |], let us consider two competing quan- 
tum Hall liquids in the presence of a long-range random 
potential, with the magnetic field determining volume 
fractions (close to 1/2), as depicted in Fig. [l](a). The 
"white" regions represent the phase with (quantized) Hall 
conductivity o~\ , such that if this phase percolates freely 
throughout the sample, the system will be on the Hall 
plateau with a xy = o\. Similarly, the "black" regions 
have local Hall conductivity o~ 2 . Near the percolation 
threshold, transport is controlled by quantum tunnel- 
ing through the various saddle points (vertices), one of 
which is shown in detail in Fig. [j](b). In the presence 
of an external electric field, the non-equlibrium Hall cur- 
rents flow inside the white and black regions. (Relation 
of this "bulk" transport description to more conventional 
Landauer-Buttiker formalism will be discussed later on). 
We assume that temperature is so low that scattering 
inside the bulk regions is absent. In this case, due to 
the continuity conditions, the currents are not allowed 
to cross the edges and must flow along them focusing at 
the saddle point Jj]]. Let us define the currents entering 
and leaving the white quadrants I and I' and those in 
the black quadrants J and J' as shown in Fig. |l|(b). For 
any distribution of the current densities across the quad- 
rants, one can relate the net currents to the four electric 
potentials on edges V\ through V4, as given by 
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As one can easily check, current conservation is met only 
if V 2 — Vi = V3 — V4. Thus we have in fact 1 = 1' and 
J = J', i.e. there is a unique white-to- white current and 
a unique black-to-black one for each saddle point. 

The pattern of the current distribution in the sample 
is determined by the "white" and "black" currents J; 
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and Ji, for each vertex i in the sample (Fig. 0(a)). This 
in general can be quite complex. Let i and j be the 
average partial current densities for the white and black 
phases respectively. The total current density is written 
j Hot = i + j- The central rigorous result in this paper is 
the following theorem: 

Theorem: If and only if i _L j, the effective conductiv- 
ity tensor a obeys the semi-circle relation in Eq. ([!]) . This 
theorem reduces the semi-circle relationship between a xx 
and a xy to a purely geometrical property of current dis- 
tributions in the two phases. 

Proof: Let us introduce, in addition to the partial 
current densities i and j, the average "white" and "black" 
electric fields e and /, related to the former by 



= cr 1 e-e, j = a 2 e-f, 



(3) 



where e is the matrix of clockwise rotation by 90°. These 
relations are valid only if scattering is absent everywhere 
except at separated points (vertices). One can obtain 
them, e.g., by drawing a vertical line across the sample 
as shown in Fig. |](a). Adding total voltage drop along 
it, and the total current across it, for white and black 
regions separately, one obtains i x — o\e y and j x = (J<xf y - 
Repeating this procedure for a horizontal line yields the 
remaining relations in Eqs. (^). Now suppose i and j are 
perpendicular to each other. This can be written as 



0e ■ f, 



(4) 



where (3 is some dimensionless parameter. The total cur- 
rent density and field are then given by jtot = i+ j, e*tot = 
e + /. Using Eqs. (||) and (||), we can relate j to t to e to t 
directly, jtot = &cfi ■ &tot > where the effective conductivity 
tensor <7 e g is given by 
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It is easy to show that Eq. (j|) implies 
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indepedent of the value of (3. This is equivalent to the 
semi-circle relation ([!]) . It is also easy to see that if the 
condition i _L j is not satisfied, we would have (3 — f3T, 
where T is some rotation matrix. In this case, the de- 
terminant in Eq. (^) would depend on (3. Thus i _L j is 
also a necessary condition for the semi-circle relation (Q) . 
Q.E.D. 

We now verify analytically that i _L j holds for an inco- 
herent (classical) system at high T. This regime is char- 
acterized by a relatively short phase-coherence length, 
l(f> ^ £pj where £ p is the percolation correlation length 
determining the typical distance between vertices. In this 
case, scattering at a vertex has a local nature, the black 



and white currents at each vertex i are uniquely related 
by 



(7) 



where parameter on depends on the scattering proper- 
ties of this vertex For our purposes, it suffices to 
consider this quantity as a phenomenological parameter 
which fluctuates from vertex to vertex. As one can prove, 
Eq. (0) and the condition that the sum of currents enter- 
ing each (black or white) region must be zero, uniquely 
define all the currents {Ii, Ji}, given the avaerage cur- 
rent density jtot- Assuming that chirality is absent in 
the neighborhood of all the saddle points, parameters Ui 
must change sign when the magnetic field is reversed 



an(B) = -Oi(-B). 
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We now show that this condition leads to i _L j. Let 
{Ii, Ji} be the current distribution for some magnetic 
field B and some jtot) with i and j being the partial 
current densities. At the opposite direction of the mag- 
netic field, B' = —B, there exists a solution J[} such 
that 1[ = Ii and J[ = — Ji. The partial current densities 
are i' — i and j' = — j. On the other hand, since the sys- 
tem is isotropic and non-chiral at large length scales, the 
pair of vectors can be obtained from the pair i,j by 
rotation and reflection. These conditions are compatible 
only if i _L j. The proof given can be easily generalized to 
the more complex case where each separate vertex has a 
random chirality, but the system as a whole is non-chiral. 

Before going on to show i _L j in the quantum trans- 
port regime, let us discuss a rather general issue of the 
relation between bulk and edge-state transport pictures 
when applied to our problem. The bulk current descrip- 
tion employed above J^]- Q differs from the conventional 
edge-state formalism which assumes that current flows 
directly on the edge. This discrepancy comes from the 
use of single-electron picture in the latter theory: non- 
equilibrium current in the bulk is impossible since both 
concentration of electrons and built-in electric field (in 
our case, that from random potential) are fixed. Adding 
electrons to the edge by shifting chemical potentials re- 
sults in this picture only in an additional edge current. 
But since real electrons do have charge, accumulation of 
electrons at edge-states (and of localized quasiparticles in 
the bulk) generates an electric field E(r) inside the bulk 
regions in Fig. 0(a), causing additional Hall current j(r) 
to flow in the bulk ||. Using Lorentz transformation, one 
can show that a homogeneous interacting electron system 
subject to a homogeneous electric field acquires a drift 
velocity equal to the classical value for a single electron. 
In the spirit of linear response, it is natural to assume 
that the same is true locally for a weakly inhomogeneous 
case: j(r) = cr xy (r)e ■ E{r), where cr xy (f) = en{r)c/B. 
Our next assumption is that quasi-particles in the bulk 
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are localized by short-range impurities, so that the effec- 
tive concentration n which contribute to the Hall current 
is given by that of the incompressible liquid. As a result, 
integrating electric field and the current density along 
any curve connecting the opposite sides of an angle in 
Fig. 0(b) gives relations of Eqs. (0). 

Although the above arguments suggest that the local 
a xy picture employed here is more appropriate for de- 
scribing current distributions in the case of long-range 
electron interaction than the conventional edge-state for- 
malism, it turns out that, in fact, both descriptions will 
yield identical results, in the absence of scattering in the 
bulk. This is because, in the edge-state formalism, the 
current density can be written as 

jedge = tp(r)e ■ Vo- X y, (9) 

where etp(r) is the chemical potential function. In the 
local a xy picture, the corresponding current density is 
given by 

Jioc = -cr xy e ■ Vip(r), (10) 

where (p(r) is now the electric potential function. For 
step-like distribution of a xy [from white region to black 
region in Fig. |l|(a)], Eq. (^|) yields a delta-function on the 
edge with an integrated current (a 2 — ci)^ e dge- As is easy 
to check, the difference between the two current densi- 
ties @ , ([l0]) is a curl which yields zero upon integration 
along any line throughout the whole sample, when the 
two cp functions in the two formalisms are identified to 
one another. Thus we are free to choose either formalism 
according to our convenience in our problem. 

We now proceed to show that i _L j for coherent quan- 
tum percolation systems. At very low temperatures, dif- 
ferent clusters of the same liquid couple to each other 
through quantum tunneling to form a coherent wavefunc- 
tion whose size is larger than the classical cluster size £ p 
and diverges when approaching the percolation thresh- 
old. Quantum coherence implies that currents at a given 
saddle point can depend on those at other ones, so that 
the local relation Eq. (^) is no longer valid. Moreover, it 
is no longer sufficient to use currents {ii, Ji} to describe 
a given state of the system, since complex quantum am- 
plitudes must be used instead. Thus we need to resort to 
a particular quantum percolation model. Our choice is 
the fully quantum coherent network model first suggested 
by Chalker and Coddington ||, which has proven quite 
successful in the study of localization length exponent in 
quantum Hall systems ||, (7). In this model, a complex 
amplitude Ak is assigned to each bond of an oriented reg- 
ular square network [see Fig. ^(c)]. At the vertex marked 
Q in Fig. ^(c), the four quantum amplitudes are related 
by a scattering matrix given by 

/A 4 \_/e t ^ \ /cosh 7 sinh 7 \ /e* 

V e^J Vsinh7 cosh 7/ \ e i ^)\A 2 ) 

(11) 



where <pk are random phase factors which model the ef- 
fect of random path lengths in an actual random geom- 
etry such as in Fig. |l](a), 7 is the tunneling parameter 
which, for the sake of simplicity, is set to be the same for 
all the vertices. For vertices oriented like the one marked 
T in Fig. |l](c), the bonds denoted Ai,...,Ai should be 
rotated by 90° with respect to those with orientation O. 
The limits 7^0 and 7 — > 00 correspond to 100% perco- 
lation of the white or the black regions, respectively. The 
value 7 C = ln(l + defines the quantum percolation 
threshold |(|. The system of Eqs. ( pi] ) written for all the 
saddle points determines the distribution of the quantum 
amplitudes, i.e. the wavefunction. 

In order to complete the construction of our formal- 
ism, we need to adapt this model to the study of trans- 
port properties. Original quantum percolation model Q 
used the edge-transport theory as a starting point, its 
classical analog being the "network model" suggested in 
Ref. ||. The squares of amplitudes \Ak\ 2 were identi- 
fied with the edge currents entering and leaving the ver- 
tices as shown in Fig. [j] (c). Unitarity of the scattering 
matrix in Eq. ([ll]) takes care of current conservation: 
l^il 2 ~ \A 2 \ 2 = \A 4 \ 2 - \A 3 \ 2 . Obviously, the standard 
edge-transport theory cannot be applied consistently to 
the quantum case since the chemical potential it deals 
with cannot be introduced for separate bonds belonging 
to the same wavefunction. Hence we identify \Ak\ 2 with 
the fictious edge current (a 2 — cri)Vk discussed above in 
a formal sense, with Vk being the electric potential of the 
bond. Then, as one can see from Eqs. (|^), the physi- 
cal currents in the white and black quadrants [Fig. |l|(b)] 
are related to the differences in the edge currents on the 
adjoining edges, given by 

I = a, (\A 2 \ 2 - {A,] 2 ) , J = a 2 (\A 3 \ 2 -\A 2 \ 2 ), (12) 

apart from a constant factor o\ — a 2 . Since we study 
small non-equilibrium currents everywhere in this work, 
Eqs. ( |l2| ) can be linearized 

/ = o 1 \Re{A 2 a' 2 ) - Re{A x a\% 

J = a 2 [Re{A 3 a* 3 ) - Re(A 2 a* 2 )}, (13) 

where ak are small perturbations of the quantum ampli- 
tudes due to an applied external voltage V, while Ak are 
the equilibrium values in the absense of V. Quantities 
Re{Aka* k ) are the non-equilibrium electric potentials. 

The remaining task is to solve the system of linear 
Equations ( pT| ) numerically with proper boundary con- 
ditions and find the average black and white current 
densities using Eqs. (fl3|). We chose standard periodic 
boundary conditions on the left and right boundaries of 
a rectangular sample in Fig. |l](c). Since Eqs. ( p"T| ) are 
linear, they can be solved separately for Ak and ak, but 
with different boundary conditions. We let all equilib- 
rium amplitudes Ak on the incoming bonds on the upper 
and lower boundaries to be unity. For the same sample, 
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we let incoming amplitudes afc on the lower boundary to 
be +V/2 and those on the upper boundary to be —V/2. 
This simulates the situation where a small two-terminal 
voltage V is applied across the Corbino disc-shaped sam- 
ple attached to two metallic contacts. Indeed, as shown 
by one of us when studying the role of contacts in the 
non- uniform Hall conductors fJ, a potential step near 
the contact is possible if and only if local a xy increases to 
the right when looking into the sample from the contact. 
This is the case for the bonds arrowed as outcoming ones 
in Fig. |l|(c) (<ti < 02 is implied). The incoming bonds, 
on the contrary, must have potentials of the contacts. 
This is completely analogous to the "ballistic" boundary 
conditions in the conventional edge-state transport the- 
ory where the incoming (leaving a contact) one-electron 
states are assumed to have the chemical potential of the 
contact. 

For each realization of {4>k}, we can extract the cur- 
rents densities i and j by adding up the corresponding 
currents at the vertices, according to Fig. |l|(c). Namely, 
i x and j y are found by averaging I and J currents, respec- 
tively, over all the vertices with orientation fl (Fig. [j](c)). 
Similarly, i y and j x are extracted from the vertices T. 
The total number of rows of vertices in either direction 
is even. The system size L x x L y ranges from 20 x 20 
to 50 x 82, and we averaged over 2,000 realizations to 
find the ensemble averaged (i) and (j) currrents. Wc 
note the importance of ensemble averaging: in the fully 
quantum mechanical case studied here, there is no sim- 
ple self-averaging even in a large sample. Self-averaging 
occurs, in a physical system, only on the length scale of 
an inelastic phase-breaking length, which at low temper- 
atures can be very large, the effect of which is replaced 
in our simulation by a finite system size. Our result for 
7 = 7c = ln(l + y/2) can be summarized as 



in the regime of quantum transport. 
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FIG. 1. (a) Schematic illustration of the random quantum 
percolation model for describing the transition between suc- 
cessive quantum Hall plateaus, (b) A single saddle point mag- 
nified, (c) The Chalker-Coddington network model (cti < az). 
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cosa = v i KJ J = 0.01 ± 0.02. 
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This shows that the white and black currents are indeed 
perpendicular to each other to a high degree of accuracy. 
We have checked that this null value for cosa persists 
when the size of the sample is changed. Neither did we 
find any dependence of cos a on the sample shape for 
any L y /L x > 1.3. Finally, we made sure that cosa, in 
such samples, is insensitive to changes of 7 around the 
critical value j c . The interval 7 in which the ratio of 
|z|/|j| changes by 2 orders of magnitude was studied. (A 
weak shape dependence at 7 = "f c was revealed in short 
samples: at L y /L x — 0.3, cosa = 0.14 ± 0.02. We in- 
terpret this as an artifact of closely placed metallic con- 
tacts; since when the localization length becomes smaller 
than the sample size as we take 7 7^ j c , cosa restores its 
null value.) In combination with the Theorem discussed 
before, we have thus successfully demonstrated that the 
universal semi-circle relation between a xx and a xy holds 
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